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Abstract 

We study the gradient flow model for the Landau-de Gennes energy functional for nematic 
liquid crystals at the nematic-isotropic transition temperature on prototype geometries. We 
study the dynamic model on a three-dimensional droplet and on a disc with Dirichlet bound¬ 
ary conditions and different types of initial conditions. In the case of a droplet with radial 
boundary conditions, a large class of physically relevant initial conditions generate dynamic 
solutions with a well-defined isotropic-nematic interface which propagates according to mean 
curvature for small times. On a disc, we make a distinction between “planar” and “non-planar” 
initial conditions and “minimal” and “non-minimal” Dirichlet boundary conditions. Planar ini¬ 
tial conditions generate solutions with an isotropic core for all times whereas non-planar initial 
conditions generate solutions which escape into the third dimension. Non-minimal boundary 
conditions generate solutions with boundary layers and these solutions can either have a largely 
ordered interior profile or an almost entirely disordered isotropic interior profile. Our examples 
suggest that whilst critical points of the Landau-de Gennes energy typically have highly local¬ 
ized disordered-ordered interfaces, the transient dynamics exhibit observable isotropic-nematic 
interfaces of potential experimental relevance. 


1 Introduction 

Nematic liquid crystals are classical examples of mesophases between conventional solid and liquid 
phases; they are anisotropic liquids with preferred directions of molecular alignment, these direc¬ 
tions being referred to as ’’directors” in the literature mm- In other words, nematic liquid crystals 
are complex liquids with orientational order. Nematics in confinement are an exciting source of chal¬ 
lenges for mathematicians and practical scientists alike. We study dynamically metastable nematic 
configurations with interfaces in prototype geometries within the Landau-de Gennes (LdG) theory 
for nematic liquid crystals. In this framework, the nematic state is described by the Q = {Qij}- 
tensor, a symmetric, traceless 3x3 matrix, interpreted as a macroscopic measure of the nematic 
order mm- The LdG energy functional comprises a bulk potential, determining nematic order as 
a function of temperature, and an elastic energy density which penalizes spatial inhomogeneities. 
We work at the nematic-isotropic transition temperature, where both the isotropic and nematic 
phases are minimizers of the bulk potential, and with the one-constant elastic energy density in 
the limit of vanishing elastic constant. This limit, describing macroscopic domains (see 0) with 
length scales much larger than the nematic correlation length, is studied in detail in the context of 
energy minimizers in [4]. 

We adopt the gradient flow model to describe the nematodynamics in the absence of fluid flow at 
a constant temperature. Gradient flows are evolution equations driven by a decreasing energy [5]. 
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Our model is derived from the LdG energy with the L^-norm as the dissipation mechanism and is 
described by a system of five coupled nonlinear parabolic partial differential equations. Standard 
theory for parabolic systems shows that we have a unique “dynamic” time-dependent solution for 
physically relevant initial and boundary conditions. Gradient flows have been used in the context 
of LG dynamics. For example, in mm the authors study a one-dimensional gradient flow model 
and the effects of biaxiality and elastic anisotropy. In [8], the authors study isotropic-nematic front 
propagation using the method of matched asymptotic expansions within the more general Beris- 
Edwards theory for nematodynamics [9]. In particular, they account for fluid flow and the coupling 
between fluid flow and nematic order. They derive evolution laws for the velocity field, the director 
field of nematic alignment and the isotropic-nematic interface but without any special attention to 
the effects of boundary conditions and initial conditions. We work in a simpler dynamical framework 
with no fluid flow but with focus on how the dynamics is affected by the choice of boundary and 
initial data. 

At the nematic-isotropic transition temperature, the LdG bulk potential bears strong resemblance 
to the Ginzburg-Landau (GL) potential in superconductivity. In our first model problem we study 
a three-dimensional droplet with Dirichlet radial boundary conditions. We use the concept of 
“normalized energy” for the GL gradient flow model in to prove that isotropic-nematic 

interfaces propagate according to mean curvature in certain model situations. The long-time dy¬ 
namics is described by an explicit critical point of the LdG energy - the radial hedgehog (RH) 
solution [El ISl E]. The RH solution has perfect radial symmetry, with perfect radial nematic 
alignment and an isolated isotropic point at the centre, referred to as a point defect in the liter¬ 
ature m- We focus on the interplay between initial conditions and transient dynamics followed 
by convergence to the static RH solution using four different representative initial conditions. The 
transient dynamics has some universal features which may have experimental repercussions. 

In Section our second model problem focuses on dynamic solutions on a disc with Dirichlet 
conditions, subject to two distinct types of initial conditions: planar and non-planar. Planar 
Q = {Qijj-tensors, = 1,2,3, have zero Qis and Q 23 components and non-planar Q-tensors do 
not. Using standard techniques, we prove that planar initial conditions evolve to planar dynamic 
solutions which have an isotropic point at the centre of the disc for all times. These solutions develop 
an isotropic-nematic interface which propagates inwards and is arrested at the origin. Non-planar 
initial conditions, including small “non-planar” perturbations of planar initial conditions, converge 
to a universal non-planar profile. The “small-time” dynamics are almost indistinguishable from the 
planar case; however, the interface collapses at the origin and the dynamic solution escapes into 
an entirely ordered non-planar state. We track the transient dynamics and numerically compute 
quantitative estimates for the “persistence time” of the interface. 

These numerical results are complemented by some analysis for radially symmetric planar critical 
points of the LdG energy that have been reported in [15] for low temperatures. We generalize some 
of the results in [15] to the nematic-isotropic transition temperature and use these critical points to 
construct radially symmetric and non-symmetric initial conditions for the numerical simulations. 
The different types of initial conditions suggest that the transient dynamics have universal features 
independent of the symmetry or uniaxiality/biaxiality of the initial condition. Namely, in all 
cases, we have a well-defined isotropic-nematic interface as a pronounced feature of the evolution 
trajectory, which we illustrate by the eigenvalue evolution of the corresponding LdG Q-tensor and 
plots of IQ 1^ as a function of time. 

The first two model problems have minimal boundary conditions which are minimizers of the LdG 
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bulk potential. In the last section, we study two-dimensional (2D) and three-dimensional (3D) 
LdG dynamic solutions on a disc with non-minimal boundary conditions. The two-dimensional 
case can be easily understood and all dynamic solutions exhibit a rapidly growing isotropic core 
with a thin boundary layer near the lateral surface. The three-dimensional solution landscape 
is richer and the transient dynamics is sensitive to the initial condition. If the initial condition 
is planar with an isotropic-nematic interface relatively close to the centre of the disc, then the 
interface propagates towards the centre, replicating the planar dynamics in Section]^ If the planar 
initial condition has an isotropic-nematic interface relatively close to the boundary of the disc, 
the interface propagates outwards, yielding an almost entirely isotropic interior and replicating the 
two-dimensional dynamics. In all cases, we have a boundary layer to match the fixed non-minimal 
Dirichlet condition. Here, the transient dynamics is sensitive to the initial interface location, a 
feature which is missing in the model problems with minimal Dirichlet boundary conditions. 

The paper is organized as follows. In Sectionwe present the gradient flow model for the Landau- 
de Gennes energy. In Section we study dynamic solutions on a droplet with Dirichlet radial 
conditions. Section follows with emphasis on planar and non-planar initial conditions on a 
disc and Section illustrates the diverse possibilities with non-minimal boundary conditions. We 
conclude in Section with future perspectives. 

2 Preliminaries 

The LdG Q-tensor order parameter is in the space of symmetric traceless 3x3 matrices, Sq = 
{Q G : Qij = Qji^ Qii = O}. A Q-tensor is said to be (i) isotropic if Q = 0, (ii) uniaxial if Q 

has a pair of degenerate non-zero eigenvalues and (hi) biaxial if Q has three distinct eigenvalues 
pile]. A uniaxial Q-tensor can be written as Q-^ = s(n(8)n — 1/3) with 5 G M and n G 5^, a unit 
vector. The scalar, 5, is an order parameter which measures the degree of orientational order. The 
vector, n, is referred to as the “director” and labels the single distinguished direction of uniaxial 
nematic alignment mm- 

We work with a simple form of the LdG energy given by 

^[Q]= / TvQp + /B(Q)dy (1) 

where 

The variable A = a(T — T*) is the re-scaled temperature, a^L^B^C > 0 are material-dependent 
constants and T* is the characteristic nematic supercooling temperature mm- Further trQ^ = 
QijQij and trQ^ = QijQjkQki foi* L j? fc = 1, 2, 3. It is well-known that all stationary points of the 
thermotropic potential, /b, are either uniaxial or isotropic [DIISIIIZ]. The re-scaled temperature 
A has three characteristic values: (i) A = 0, below which the isotropic phase Q = 0 loses stability, 
(ii) the nematic-isotropic transition temperature, A = BA jTJC^ at which is minimized by the 
isotropic phase and a continuum of uniaxial states with s — — BjSC and n arbitrary, and (iii) 

the nematic supercooling temperature, A = S^/24(7, above which the ordered nematic equilibria 
do not exist. 

Throughout this paper we work at the nematic-isotropic transition temperature, investigating 
the propagation of fronts separating the isotropic phase from the ordered equilibria in the limit. 
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LCjR^B^ 0^, where i? is a characteristic length scale of the domain Vt. We refer to this as the 
vanishing elastic constant limit for fixed values of by analogy with the terminology in | 3 ]. 

Continuum formulations are typically valid in this limit mm- We work with the gradient flow 
model associated with the LdG energy Hi and the dynamic equations are given by: 

IQI')-CIQI'Q, (3) 

where 7 is a positive rotational viscosity, QQ = QijQjk with i^j^k = 1,2,3 and I is the 3x3 identity 
matrix. The system comprises five coupled nonlinear parabolic partial differential equations. 
We recall a basic result about the existence and uniqueness of solutions for such gradient flow 
systems: 

Proposition 1. Let Lt C be a hounded domain with smooth boundary, dLt. Given a smooth 
fixed boundary condition Q(x,t) = on dLl and smooth initial condition, Q(x,0) = Qo(x), 

the parabolic system has a unique solution, Q(x,t) G C^{Ll) for all t > 0. 

Rroof. The existence of a solution is standard; see m for a proof. From mm , we have the dynamic 
solution is bounded for all times with |Q(i‘,t)| < BjSC for t > 0. The uniqueness result 

follows from an immediate application of Gronwall’s inequality to the difference = Qi — Q 2 of 
two solutions, Qi and Q 2 , subject to the fixed boundary condition and the same initial condition. 
In particular, Qci(x, t) = 0 on dQ and Qci(x, 0) = 0 for x G Lt. One can then show that Qd(x, t) = 0 
for X G and for all t > 0. □ 


7 Qt = LAQ -AQ + B 
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3 Front propagation on three-dimensional spherical droplets 


Our first example concerns nematic droplets. Let Q be the unit ball in three dimensions; Q := 
|xgM^;|x|<1}. We work with a uniaxial Dirichlet boundary condition 


^ . I 

'*'■=30 '■*■■“3 


(4) 


where f is the 3D radial unit vector. is a minimizer of the bulk potential fs in Q. For 
illustration, we first work with uniaxial radial initial conditions that have a front structure such as 


Q(x,0) 


0 0 < |x| < ro 

Qb ro < |x| < 1, 


( 5 ) 


for some ^ < tq < 1. We refer to these as “radial hedgehog” type initial conditions by analogy with 
the static radial hedgehog solution as described in Section [El [m E]. We are interested in the 
qualitative properties of dynamic solutions of Q subject to these initial and boundary conditions, 
such as front propagation and transient dynamics. Looking for “dynamic” radial hedgehog type 
solutions, we work with an ansatz of the form 


Q(r,f) = h{r,t) 



( 6 ) 
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where h : [ 0 , 1 ] x [ 0 , oc) ^ M is the scalar order parameter that only depends on r, the radial 
distance from the origin, and time. On substitution into Q, we have a solution of the form Q if 
the scalar order parameter /i is a solution of 


ryhf = hrr + -hr - ^ (/i+ - h) (2h - /i+), 


( 7 ) 


with ^ and L = The boundary conditions are /i(0,t) = 0, 

the initial condition is 


/i(r, 0) 


0 0 < r < ro 

ro < r < 1. 


for alH > 0 and 


( 8 ) 


The evolution equation for /i in Q is simply the gradient flow model associated with the energy 
functional 


m 



{h-hlY\ 

j 


dr. 


( 9 ) 


Given a smooth solution h{r^t) G (7^([0,1] x [0,oo)), we can appeal to the uniqueness result in 
Proposition to deduce that Q is the physically relevant solution and hence, the five-dimensional 
evolution problem in Q reduces to a single evolution equation for a scalar order parameter. 

In [To], the authors study a closely related problem for front propagation on three-dimensional 
balls in the Ginzburg-Landau framework. They rigorously prove that for suitably defined initial 
conditions (as in Q) with appropriately bounded energy, the front propagates according to mean 
curvature. Our governing equation Q is similar to that studied in Section 3 of m, however we 
have an extra term: —6h/r‘^ in Q. In particular, we cannot quote results from uni and m without 
verifying that the key inequalities are unchanged by the additional term for L sufficiently small. In 
the next paragraphs, we verify the necessary details to reach the desired conclusion. Let p{t) be 
the solution to 


^ = --, p(0) = roG (1/2,1), 

at p 


( 10 ) 


or alternatively, p{t) = y/rg — 4L We define Ti = | (tq — |). This is the first time for which 
p{t) = 1/2 and is independent of L. Next, let 


/(g t) 


0 0 < r < p{t) 

h^ p{t) < r < 1. 


( 11 ) 


Our goal is to show that the solution h{r^t) of 0 , subject to suitably defined initial conditions, 
resembles the function /(r, t) for T < Ti, in the sense that 


[ \h{r,t) - /(r,f)| dr 0 

Jo 


as L —> 0 . 


As in sni, the key step is to define a weighted energy as shown below: 


( 12 ) 




+ 


w 


(/l+ 


w) 


Vl 


di?. 
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where 


(13) 


w{R^ t) = h{R + p(t)R)^ ~p{^) ^ R ^ ^ — p(^)5 '7" > 0, 

and 0 (i?, r) is a weight function 

2 


r) = exp 


P . 


1+5 


(14) 


In particular, w {R, r) is a solution of 


\/l/ . Qw^/L 3 1 \ n 

V Lwr -— {(PWR)^ + ———^ + —u;(/i+ - w){2w - h+) = 0 . 

(p " [R + pY Vl 

We follow the steps in Proposition 3.2 of m to show that 


(15) 


^E4 w]{t) < P P{R,r)wl dR+^L^/^hl, 

UT O 'J J _p(^\ p 


(16) 


for T < Ti. By contrast, in sni, the weighted energy in the GL-framework is strictly decreasing. We 
have lesser control on the weighted energy but for L sufficiently small, the inequality (16) suffices 
for our purposes. 

Next, we define an interface energy which yields lower bounds for the weighted energy: 


2 ^ 

g{s) = —j^ w{h+- w) dw, 


(17) 


hi . 


and so g{hj^) — is an interface energy associated with an isotropic-nematic front. Further, let 

0 ~p('i) < i? < 0 


viR) = , 

0 <i?<l — p{t). 

We can adapt a lemma in m to show that: 


(18) 


Proposition 2. If for some smooth function w, 


f 


\g(w) - g(v)l ds < and E^[w] < Ci, 


(19) 


where 0 < cr < 1/4 and a = p{Ti)/2y^, then E(f)[w]{r) > g{h-^) — C 2 L^^‘^ ^ — CsL^^ for r < Ti 
and positive constants Ci,C 2 ,C 3 independent of L. 

It remains to construct initial conditions w{R^0)^ which satisfy the hypothesis of Proposition]^ 
The construction is parallel to that in Equation ( 1 . 22 ) of [ 2 Tj and we give a statement for com¬ 
pleteness. 
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Proposition 3. Define the funetion 


a{R) = 


h. 


1 + exp[—^/3h^R] 
At T = 0, R = r — ro and for L sufficiently small, define 
(h^ 


w (R) — < 


(/i+ - a (^)) {R - 2LV4) /lV4 + h+ 

(a) 

-^) (i? + 2Li/4) /lV4 


a 


a - 


R > 2Li/^ 

lV4 <R< 2LV4 
-LV4 <R< lV4 
-21^/^ <R< -L^/^ 

i? < -2LV4. 


Then E^[w] < g{h^) + CLf!'^ for a positive constant C independent of L, os L —>■ 0. 
Finally, we have by analogy with the main theorem in m- 


( 20 ) 


( 21 ) 


Proposition 4. Let 0 < a < ^ and assume that E(p[w]{0) < g{h^) + ciL‘^^ for some eonstant 
Cl > 0 independent of L and that for a as above, 

j‘'y{w{R,t)))-g{v)\ di?<^L“ (22) 

Let Te he the first time for whieh 

P \g{w{R,T,)) - g{w{R,0))\ di?= (23) 

Then > min(Ti,C) for some positive eonstant C independent of L as L ^ 0'^. In other words, 
we have 

J\{wiR,T))-g{v)\ dR<^L- (24) 

for all T <T^ and is of order one. 

The proof follows verbatim from HU. Equipped with a weighted energy, estimates for the rate 
of change of the weighted energy and bounds for the weighted energy along with suitable initial 
conditions, we adapt arguments from Theorem 3.1 of m to prove: 


Theorem 1. Let D be the unit ball in Let hp {r,t) denote the solution of the evolution equation 
0 suhjeet to the fixed boundary eonditions and an initial eondition with an interfaee strueture and 
appropriately bounded weighted energy: 


f 


i’oir) 


2w^ 




+ 


w‘^ (/i+ — wf 


' dr < g{h+) + 


(25) 
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with r independent of L and fjQ (r) = ^ exp 
been defined in Proposition\^ we have 


-2 . Then for any T < T^, where 


has 


lim 


’ r 

/ \hi{r,t) - f{r,t)\r^drdt = 0 

Jn 


(26) 


where f{r^t) has been defined in (11). 


Comments on the proof: The condition (25) is equivalent to the bound in Proposition]^ which 
can be realized by initial conditions with “efficient interfaces”. The key ingredients are the rate 


of change of the weighted energy in (16), the lower bound for the weighted energy in terms of 


g{h^) in Proposition and the upper bound in (25). These ensure that the system does not have 
sufficient energy to create additional interfaces away from R = p(t) and h is effectively constant 
(either /i = 0 or /i = ) away from R = p(t). Given hp (r, t), we obtain the unique solution 

Ql(x,t) of the evolution equation (|^, subject to the boundary and initial conditions, given by 

Ql (x,i) = /iL(|x|,t) (f ®f- I). 


3.1 Numerical simulations on the sphere 

We numerically compute solutions of the gradient flow system © on a three-dimensional droplet, 
with the fixed boundary condition in (|^ and various types of initial conditions. From the nu¬ 
merical results in m, the radial hedgehog (RH) solution is the global minimizer of the LdG energy 
for this model problem; the RH solution is a uniaxial solution of the form H = 5(r) (f 0 f — |) 
where 5(0) = 0 and 5(r) > 0 for r > 0. In particular, s rapidly interpolates between 5 = 0 and 
the boundary value of 5+ = BjZC over a distance proportional to the nematic correlation length, 
oc ^JLCIB‘^^ and the localized region of reduced order near r = 0 is referred to as the ’’defect 
core”. We expect the long-time dynamics to converge to the RH solution for all choices of initial 
conditions. However, we are equally interested in the transient dynamics and the dynamic persis¬ 
tence of isotropic-nematic interfaces. In what follows, we look at four different initial conditions: 
uniaxial initial conditions within the remit of Theorem uniaxial initial conditions outside the 
scope of Theorem]^ biaxial initial conditions and initial conditions that break the radial symmetry 
of the order parameter. 

Let R be the radius of the droplet and we non-dimensionalize the system (j^ by setting i = x = 









^ to yield 


^ ( ^Qii + 2 C(( 5 ii + Q22 + Qi2 + Q11Q22 + Qi3 + ^23)^11 




+ Qi 2 + Qi 3 ~ 2Q22 ~ ‘^ QiiQ 22 — 2(523) ) 5 


(27) 


— A (522 — ^ ( ^(522 + 2 ( 7 (( 5 ii + Q22 + Qi2 + Q11Q22 + ( 5 i 3 + ^23)^22 


B 


~^iQl 2 + Q22 + Q23 ~ 2 ( 5 ii — ‘2Q11Q22 — 2 ( 5 i 3 , 


(28) 


— A(5i2 — ^ ( 2iQi2 + 2C{Q‘ii + (522 + Qi2 + Q 11 Q 22 + (5l3 + ^23)^12 


-S((5 ii (5 i 2 + Q12Q22 + Q13Q23) ) , 


(29) 


— A( 5 i 3 — £ ( ^( 5 i 3 + 2 ( 7 (( 5 ii + Q22 + Qi2 + Q11Q22 + ( 5 i 3 + ^23)^ 


13 


— B{Qi 2 Q 23 — Q22Q13) ) 5 


(30) 


—— AQ23 — ( ^Q 23 + 2 ( 7 (( 5 ii + Q22 + Qi2 + Q11Q22 + ( 5 i 3 + ^23)^23 


—B{Qi 2 Qi 3 — Q23Q11, 


(31) 


where A denotes the Laplacian with respect to the re-scaled coordinate x. (In what follows, 


we drop the bars from the dimensionless variables). We take = 10 L = 


B — 0.64 X WrNjvfr^ C — 0.35 x IQ'^N/rrr and A = throughout the paper and work with 
either L = 0.05 or 0.01 [l3j. 

Here, and in the subsequent sections of this paper, the system of reaction-diffusion equations 


(27 - 31) is solved as follows. The unit ball is embedded into the unit cube [—1,1]'^ which is 


discretised with a uniform cartesian grid with spatial resolution h. We implement a special case 
of an immersed boundary method (see for example, [22]), and apply the boundary conditions at 
all discrete points within distance h/2 of the boundary. For interior points, the solution satisfies 
the system ( 27]|^ ) and in the exterior of the physical domain, we solve the simple heat equation 
(i.e. take A = B = C = 0) and use periodic boundary conditions on the cube. This setup makes it 
simple and efficient to use higher order and spectral schemes for spatial derivatives. Timestepping 
is accomplished with a standard fourth-order Runge-Kutta scheme. Simple finite difference schemes 
are also implemented to verify the results. 

The first two initial conditions are uniaxial RH type initial conditions of the form 


Q(r,0) = h{r,0) 


r- 


(32) 


Case I prescribes an initial condition, /i(r, 0), with an interface structure given by h{r) = ^1 -h tanh ^(r — tq) / 

and Case II describes an initial condition without an interface structure with h{r) = h^r. Case 
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Figure 1: |Q(r,t)p on the cross section of the unit ball at 0 = 0 for Case I and Case II, at t = 0, 
t — 0.001 and t — 0.125. The spatial resolution is h — 



Figure 2: Interface position for Case I and Case II with ro = 0.5 (black), and predicted position 
according to motion by mean curvature (blue) for L — 0.05 (dashed) and L = 0.01 (solid). In Case 
II an interface quickly develops so the two curves are indistinguishable. The spatial resolutions for 
L — 0.05 and L — 0.01 are h — and ^ respectively. 


I is within the remit of Theorem [T] and the numerics demonstrate that the solution retains the 
isotropic-nematic interface for all times and that the interface propagates towards the origin ac¬ 
cording to mean curvature for small times, equilibriating near the origin for long times. For long 
times, the radius of the isotropic core scales, as expected, with ^/L and arises out of the saddle 
structure of Q at the origin. The dynamic solution in Case II very quickly develops an inwards- 
propagating interface separating the isotropic core at the centre from the ordered nematic state 
and then follows the same evolution path as Case I (see Figure]^. The long-time behaviour of the 
dynamic solutions for Cases I and II are indistinguishable within numerical resolution, as expected. 



Figure 3: Eigenvalues of Q(r,t), as a function of r, for Case I and Case II (dotted), at t = 0, 
t = 0.001 and t = 0.125. 
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Figure 4: Eigenvalues of Q(r, t), 
and t = 0.085. 


as a function of r for initial condition (33), at t = 0, t — 1.5 x 10 


-4 



Figure 5: Qii(r, t) — /i+((cos (/)sin0)^ ~ V^) on the cross section of the unit ball for initial condition 
(33), at t = 0, t = 1.5 X 10“^ and t = 0.085. The spatial resolution is /i = 


For Case III, we use a biaxial initial condition given by 


Q{r,0) = h{r) 


(g) f 



+ s{r) {m ® m 




(33) 


where m — (cos 0 cos 0, cos 0 sin 0, — sin0) and p — (—sin 0, cos 0,0). The function h{r) has an 
interface structure, as in Case I, and s{r) = r(l —r). This case is outside the scope of Theoremand 
we are not guaranteed the uniaxial radial symmetry of the dynamic solution. The numerics show 
that the dynamic solution quickly becomes uniaxial within numerical resolution, as demonstrated 
by the evolution of the eigenvalues of Q(r, t) in Figure]^ The dynamic solution exhibits an inwards- 
propagating interface separating the isotropic core at r = 0 from the ordered nematic state and the 
interface equilibrates near the origin. We also numerically compute the differences 

Qij{r,t) Sij\ 

IQI V |rP 3 ; 

as a function of time and our numerics show that this difference vanishes everywhere away from 
the origin, since |Q(0,t)| = 0 for large times; see Figure]^ 


Case IV breaks the radial symmetry of the initial order parameter by employing an uniaxial initial 
condition of the form (32) with 


, , ^ . B ( , f sin^ 6 >(cos^ 0 + 4 sin^ 0 ) + 2r^ cos^ 0 — 0.5 

h{r, 0, 0 ) = ^ 1 + tanh -^^ 


6C 




(34) 


Here, the initial interface is ellipsoidal in shape. The interface becomes circular and the subsequent 
dynamics are indistinguishable from Case I as seen in Figure 
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Figure 6 : |Q(r,t)p on the cross section of the unit ball for initial condition (34), at t = 0, t = 0.05 
and t = 0.1. The spatial resolution is h = ^. 


These examples illustrate that whilst the static RH solution has a localized defect core of reduced 
order near the origin that may not be experimentally observable, the transient solutions exhibit 
well-defined isotropic-nematic interfaces (see Figure [^. These interfaces propagate towards the 
droplet centre and may be experimentally observable. 

4 Fronts on a disc 

4.1 Analysis on a disc 

We take our computational domain to be the unit disc defined by 

n = {(r,0) G 0 < r < 1, 0 < 0 < 27r} , (35) 

with fixed boundary condition Q = Qc = ^ (ni (g) ni — |) on r = 1 for ni = (cos 0, sin 0 , 0 ), 
the radial unit vector. This boundary condition is purely uniaxial and is a minimum of the bulk 
potential. We study dynamic solutions of the parabolic system, (§, subject to Q = Qc on r = 1 
with different types of initial conditions. We note that such solutions also survive as translationally 
invariant solutions, independent of z, on a cylinder with free boundary conditions on the top and 
bottom surfaces. 

We first present some heuristics based on critical points of the LdG energy on a disc subject to 
this Dirichlet condition. The gradient flow model dictates that dynamic solutions evolve along a 
path of decreasing energy, converging to a critical point of the LdG energy [5]. Hence, the long¬ 
time behaviour can be predicted by an analysis of the corresponding stationary problem. In [ 3 ], 
the authors present a general analysis of LdG energy minimizers on 3D nice domains (see in 
the L ^ 0 limit. Based on their analysis, the minimizers converge strongly in (D; Sq) to a 
limiting harmonic map of the form Q = s (n* (g) n* — |) such that 5 = 0 or s = B/3C a.e. (so 
that Q is a minimum of fs) and the director n* is a solution of the harmonic map equations 
An* -h |Vn*pn* = 0 (also see [23] for recent work on planar domains). The convergence is shown 
to be uniform away from the singularities of the limiting harmonic map, which need not be unique. 
There are at least two harmonic maps on a disc with the boundary condition n* = (cos 0, sin 0, 0) 
on r = 1 [ 23 ] , 


ni = (cos 0 , sin 0 , 0 ) , and 112 


2x 2y 1 — 

1 -h ’ 1 -h ’ 1 -h y 


(36) 


We conjecture that there are at least two competing limiting harmonic maps for this two-dimensional 
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problem, defined in terms of ni and 112 above: 


(37) 


Qi = s 


^ni (g) ni - 



and Q 2 


B 


^n2 (8) n2 - 



As ni is not defined at r = 0, Qi must have an isotropic point at r = 0, with s B/3C rapidly 
away from r = 0. We refer to Qi as being the two-dimensional planar radial hedgehog profile. 
However, n 2 has no singularity on and Q 2 does not have an isotropic core. We predict that 
the dynamic solutions of Q, subject to Q = Qc on r = 1, converge to Qi everywhere away 
from the origin (where ni is singular) if escape into the ^-direction is not allowed. Equivalently, 
|QL(r,t) - Qi(r,t)| ^ 0 as L ^ 0 for |r| > cr(L), cr{L) 0 where L ^ 0 (or as ^^^2 ^ 0) 
but Qi is not a stationary solution of Q. The dynamic solutions converge to Q 2 if escape into 
third dimension is allowed. Next, we have a lemma which demonstrates that escape into third 
dimension is not allowed for certain initial conditions. We refer to a Q-tensor as being “planar” if 
the components Qi 3 ^Q 23 — 0 are identically zero on and “non-planar” if not. In particular, Qc, 
is a planar Q-tensor. 


Lemma 1. Let Q(r, t) be a solution of ^ on Q, with fixed boundary eondition Q = Qc on r = 1 and 
a planar initial eondition, Q(r,0) sueh that |Q(r,0)| < (^) for r ^ Q. Then Q 13 = Q 23 = 0 

for all t > 0 . 

Proof. The proof is an immediate application of Gronwall’s inequality m- From 120111 , we have 
the following L°°-bound for the dynamic solution, |Q (r, t) | < for f > 0. The two governing 

PDEs for Qi 3 and Q 23 can be written in the form 


dQi3 

dt 

dQ23 

dt 


LAQis = T'(Q)(2 i 3 + BQ12Q23 
LAQ13 — G{Q,)Q 23 + BQ12Q13 


(38) 


where F and G are bounded functions by virtue of the bound above. We integrate by parts, 
use the fact that Q 13 — Q 23 = 0 on r = 1 and apply Gronwall’s inequality to obtain 



13 


+ <223 dE < exp 



Ql^ + dE 


= 0 , 

t=0 


so that Qi 3 = Q 23 = 0 for all t > 0 . 


(39) 

□ 


4.1.1 Radially symmetric static solutions 

We consider a particular class of planar critical points of the LdG energy on a disc, introduced 
in [15] for low temperatures (described by A < 0 in Q), referred to as radially symmetric 
solutions. The theoretical results in this section are a generalization of the results in m to the 
nematic-isotropic transition temperature. We work in a different temperature regime where the 
LdG bulk potential has two equal energy minima, and hence we cannot a priori assume that the 
results in dSlES] apply to our case. These solutions are labelled by two order parameters u and v 
that only depend on the radial distance from the origin. We can perform some explicit analysis for 
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these critical points which are good examples of planar initial conditions. We first show that there 
exist radially symmetric planar {u^ 'u)-critical points of the LdG energy on a disc at the nematic- 
isotropic transition temperature, which are unstable in the sense that the second variation of the 
LdG energy is negative for admissible perturbations. We use these critical points to construct 
planar initial conditions for the LdG gradient flow model. We then use numerical simulations 
to corroborate our heuristics that the corresponding dynamic solutions quickly develop a radially 
symmetric isotropic-nematic interface (even when the initial data is not radially symmetric) and 
converge to Qi everywhere away from the origin. Next we investigate the effect of small “non- 
planar” perturbations of the planar initial conditions and show that the solution follows the planar 
dynamics for small but noticeable times, followed by an abrupt escape into the third dimension at 
r = 0 and long-time convergence to Q 2 . As in [15], we study the LdG Euler-Lagrange equations. 


d 2 

LAQ = —Q 


21C 


s QQ 


IQI 


and look for solutions of the form 

u(r) 


Q = —(ni (g) ni — m (g) m) + v(r) ( p (8> p — - 


I +C’IQPQ 


(40) 


(41) 


where m = (—sin 0, cos 0, 0) and p = (0,0,1). The critical points (41) are referred to as fc-radially 
symmetric solutions in [151 [26] with k = 2; we note that the authors normalize the tensors in (41) 
whilst we choose not to normalize the tensors for ease of presentation. It is straightforward to 
verify that solutions of the form (41) exist if the functions u and v satisfy the following system of 
coupled second-order ordinary differential equations 


u"(r) + 

^"(r) + 


u^(r) 4u(r) 


r 

v\r) 


u 

L 


2 ry f 2 ?;^ 

^ + ^Bv + cl- + — 


V 

L 


Bv ^ f 2'L’^ 

___l_ I _|__ 

27C 3 V 2 3 




(42) 

(43) 


with ?i(0) = 
condition Q 
(P-dH) by 


'L’'(0) = 0 with u{R) — BjZC and v{^ — —B/QC^ consistent with the boundary 
= Qc on r = 1. As in [15], we can prove the existence of a solution pair {u^v) of 
appealing to a variational problem. Define the energy 



('^ 0 ^ + 


B^ 

^ 54CL 



16 ^ ^ Vy 


B f2v^ 
~ 3L^ \~9~ 



(44) 


This is the LdG energy of the (?/,'L’)-ansatz in (41), defined on the admissible set S = : 

[0,1] ^ ^ A^(0,1), u{l) = B/3C^ 'L’(I) = —B/6C}. The proof of the 

following lemma is standard and omitted here for brevity (see US]). 


Lemma 2. For each L > 0, there exists a global minimiser {u^v) G [(7^(0,1) n(7([0,l])] x 
[(7^(0,1) n G^([0,1])] of the energy (44) on S, which satisfies the ODEs for u and v in 
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Next, we look at some qualitative properties of the (?i, 'L’)-solutions. Similar questions have been 
addressed in the recent paper [26], in the temperature regime A < 0 , with the exception of the 
monotonicity argument in Lemmabelow and we reproduce all necessary details for completeness. 


Lemma 3. Let (u,v) be a global minimizer of the energy £ in (44) subjeet to u{l) = ^ and 
'?;(1) — Then 0 < ^ and — ^ < 0 for 0 < r < 1. 

Proof. We can prove the non-negativity of u by following the arguments in [151ES] (the authors use 
the symmetry E [u^ v] = E [—u^ v] and the strong maximum principle for (42) to deduce that u > 0 
since u{l) > 0). We assume 'L’(ri) = 'L’(r 2 ) = 0 and v > 0 for ri < r < r 2 . Define the perturbation 


v{r) = 


A direct computation shows that 

E[u,v] - E[u,v] = f + 


p(r) 

0 < r < n 

r 

n < r < r 2 

[r(r) 

r 2 < r < 1. 

r Cv‘^ 

f bV 


(45) 


L 9 


3C. 


u 


(46) 


contradicting the global minimality of the pair {u,v). Next, let us assume for a contradiction that 
the minimum value of v : [ 0 , 1 ] —>■ M, say Umin < — ^ at r = tq. At r = tq, the left-hand side of 
(43) is non-negative by definition of a minimum. Ftom the maximum principle (see m, s) , we 
have 

2v‘^ 


u 


Q|^ = ^ + _< 




2 3 

< —then u^(ro) < Then for u, 

B o 


902 

B'^ 

V? 


3 

< — we have 


9C2 ) 


(47) 


2C 


B^ 


— — v“ ■ H- v'^ ■ <C - 

^min ' Q ^min ^ 54^7^ 


{B + 2Cv„ 


< 


3 

B^ 


and 


(48) 


4 ' 54(72 

SO that the right-hand side of ( [43| ) is negative, yielding a contradiction. The result for u follows in 
the same manner, using the bound for v proven above. □ 

We next show that (i/, v) are monotone functions, borrowing an idea from m- We make the 
elementary observation that > 0 for 0 < r < a, since u attains its minimum value at r = 0 . 


Ue{r) = u{r) -h e(a(r), v^{r) = v{r) + ef3{r) 

with (a(l)=/3(l) = 0. In this case the second variation is given by 


6^E 


[a,/3]:= [ 

Jo 


a2 /32 a2 


r 

+ L 


27C 3 


+ /; 


/ 

f 


C 


2 o2 I ^ 'A ^ ^ 'A 'A ^ 

—u p + —va -b —V a + 


B 


C 


. 2„2 


6 


6 


6 


2BvP‘^ 2C 2 
Bual3 ^ 2Cuvafd 


+ /; 


, 3C 2 2 
10^“ + “ 


dr. 


(49) 


dr 

dr -h 


(50) 
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In particular, 6‘^E [cr,/?] > 0 for all admissible cr,/? by the global minimality of (u^v). 


Lemma 4. Let (?i, v) be a global minimizer of the energy in (44)- Then u' > 0 for r > 0 and v' < 0 
for r > 0. 

Proof. We assume for a contradiction that u and v are not monotone, so that there exist points 
^1, ^2, ^3, ^4 with ri, r 2 G (0,1), ra G [0,1) such that 


= u\r 2 ) =0 < 0 for ri < r < r 2 

'^^(^ 3 ) = '^^(^ 4 ) = 0 'u' > 0 for ra < r < r^. 


(51) 


We differentiate the ODEs for u and v in 
over r G ['ri,r 2 ] and r G [rs^r^] to get the following equalities: 

rr2 


rr 

J ri 


r, //.o 

if” >+ 


5,».9 2 , r 

— {u Y - kUU + — 

4r L 


-(43), multiply by ru' and rv' respectively, integrate 

SC 




uu'v' 


+^(n')V +g 


{B + 2Cv) 


dr = 0, 


(52) 


and 


r . >\2 

L 5'* > ^ 

(^/)2 ^ 

- 

/\2 / 

l\2 2/ 

/\ 2 

3r L 

8lP'^ 

0 

’) + Y 

0 


r 

+ L 


uu'v' 




{B + 2Cv) 


dr = 0. 


We define the perturbations (a,/3 as follows: 
a{r) = 


0 if I/' > 0 
u' if u' < 0, 


/3(r) = 


0 if n' < 0 
v' if v' > 0. 


(53) 


(54) 


These perturbations satisfy (a(l) = /3(1) = 0, since u and v attain their maximum and minimum 
values respectively, on r = 1. Substituting (a, /3) in (50) and using (52)-(53), we obtain [a, /3] < 
0 and the required contradiction. □ 

Finally, we demonstrate that this class of critical points is unstable in the static sense, at the 
nematic-isotropic transition temperature. We consider a perturbation about the critical point in 
(H^, w,, = Q, j +eVij with \'ij = 0 on r = 1 and compute the second variation of the LdG energy 
about this critical point as shown below 

+ + + (55) 

We simply use the perturbation 

1 nOr*^ r 1 — 

(56) 


100r^(l —r^)^ 

V = (ni (g) p -I- p ® ni) 


(1 + 100r2) 


in (55). This integral is evaluated numerically using numerical solutions to (42) - (43) and 5^1 < 0 
for log^o < —1.6, as illustrated in the graph below. 
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4.2 Numerical simulations 


We numerically solve the LdG gradient flow system in ( 27pT ) on a disc subject to Q = Qc on 
r — 1. We take R to be the disc radius in the definition of the dimensionless variables and the 
parameter values (i?, L, S, C) are as in Section]^ We take L — 0.01 in this section, unless stated 
otherwise. The focus is on the distinction between planar and non-planar initial conditions. Based 


on the heuristics presented in Section [4.1.H we argue that all dynamic solutions develop an interface 
separating an isotropic core at r = 0, from an ordered uniaxial nematic phase, away from r = 0. 
For planar initial conditions, the interface persists for all times and the dynamic solution has an 
isotropic core at r = 0 whereas for non-planar initial conditions, the interface collapses at r = 0 


and the dynamic solution relaxes into the uniaxial state Q 2 in (37) for long times. 


Firstly, we consider planar (i^,'?;)-type initial conditions of the form 


Q(r.0) = ?^{n. 


mi - m (g) m) + v{r, 0) ( p (g) p - - 


(57) 


where m = (— sin0,cos0,0), p = (0,0,1) and I = ni(g)ni+m(g)m + p(g)p . Recall our boundary 
conditions enforce u = B/3C and v = —B/QC on r = 1. We let u{r, 0) and v{r^ 0) have an interface 
structure as shown below: 


«(r, 0) = y j , 0) = -^ y tanh j j . 

for various values of uq and vq. It is worth noting that the corresponding dynamic solution is given 
by 

u(r t) / I\ 

Q(r, 0) = —^ (ni (g) ni - m (g) m) + v{r, t)(p(g)p--j (58) 


if the dynamic order parameters satisfy 
JUt = u 

IVt 


. . u'(r) 4:u(r) u f B‘^ 2 ^ ^ f 2 '?;^ 


r2 I \27C 3 


nr \ V ( B'^ Bv ^ 1^2 

= v''{r) + - — + C( — + — 11- -^Bu^. 


r L V27C' 3 


2 3 


4L 


(59) 

(60) 


From Proposition!^ this is the unique solution for this model problem. In fact, we can go further 
and exploit the methods in [28] to compare the isotropic-nematic interface motion in (58)-(60) with 
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Figure 8: |Q(r, on the cross section of the cylinder at z = 0 for uq = 0.6, vq = 0.4 for the initial 
condition ( [57| ), at t = 0, t = 10“^, and t = 0.25. The spatial resolution is /i = ^ 





Figure 9: Eigenvalues of Q(r, t) for initial condition (57) with uq 
and t = 0.25. 


0.6, vq = 0.4, at t = 0, t = 10 


mean curvature motion. As in Section we cannot quote results from [28] since the dynamic 
equations (58)-(60) differ from the Ginzburg-Landau model in [28] by the additional term —Aujjf 
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for 


in (59) above. However, for L sufficiently small, this term may be controllable and in Figure 
we plot the numerically computed interface location (plot r*(t) such that |Q(r,t)p < 
r < r*{t)) and find good agreement with mean curvature propagation for small times. 

A typical solution with an initial condition of the form (57) is shown in Figure i If Uq ^ Vq^ 
then Q(r,0) is necessarily biaxial but it is hard to see the biaxial character of the initial data by 
looking at |Qp. In order to see the rapid relaxation to uniaxiality we plot the eigenvalues of the 
dynamic solution, Q(r,t), as a function of time (see Figure [^. Varying the values of uq and vq 
does not change the qualitative dynamics: Q(r,t) quickly becomes uniaxial for all choices of uq 
and vq^ within numerical resolution. The dynamic solution develops a radially symmetric interface 
separating the isotropic core, Q = 0 at r = 0, from an ordered uniaxial nematic state (away from 
r = 0) and the interface equilibrates near r = 0 for long times. We have numerically computed the 
tensor-difference, Q(r,t) - Qi, as a function of time (where Qi is defined in ([ST])) and find that 
Q(r,t) ^ Qi, everywhere away from r = 0, as expected. In particular, Qi is not a solution of the 
LdG Euler-Lagrange equations and only describes the leading order behaviour away from r = 0. 
One component of this tensor difference is plotted in Figure 


Next, we consider an initial condition of the form (57) with 


/ . B f , / r- 0.6(1+ 0.25sin5(9) \\ 

u(r, 0.0) = y tanh -L_-i j j 


12C 


ri 


( 61 ) 
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Figure 10: Qn{r^t) — — 1/3) on the cross section of the cylinder at z = 0 for the initial 

condition ( [57| ) with uq = 0.6, vq = 0.4, at t = 0, t = 10“^ and t — 0.25. 
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Figure 11: |Q(r,t)p 
resolution \s h — . 


on z = 0 for initial condition (61) at t = 0 , t = 


0.06 and t = 0.2. The spatial 


This is again a planar initial condition with an interesting star-shaped isotropic-nematic interface 
that relaxes into a radially symmetric isotropic-nematic interface. The subsequent dynamics is then 
indistinguishable from the case discussed above; this is illustrated by Figures pT] and 
Finally, we consider a small perturbation to the planar {u^ 'L’)-initial conditions, that render non- 
planar initial conditions. Let 

r, 0 ) = u{r^ 0 ) (g) n — 

where, 

n = (a /(1 — e^(l — rY) cos 0 , ^(1 — e^(l — r)^) sin 0 , e(l — r)), 

I 2 = ni (g) ni -h m (g) m, 

and p is as before. The functions ?i(r, 0) and '?;(r, 0) are as defined previously. As before, the 


I2 


+ ?;(r, 0 ) ( p 0 p - 


(62) 



Figure 12: (5ii(r,t) — — 1/3) on the cross section of the cylinder at z = 0 for the initial 

condition ( [M| ) at t = 0 , t = 0.06 and t — 0 . 2 . 
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Figure 13: Spatial dependence of eigenvalues of Q(r,t) at t = 0, t = 0.001 and t — 0.25 for initial 


condition (62) with — 0.6, — 0.4. 



Figure 14: Time evolution of eigenvalues of Q(r,t) at the origin for initial condition (62) with 
I/O = 0.6, — 0.4. 


dynamic solution quickly becomes uniaxial (within numerical resolution) irrespective of i/q '^o 
and develops a well-defined interface separating an interior region with Q = 0 near r = 0, from an 
ordered uniaxial nematic state elsewhere. This interface propagates inward but instead of being 
arrested at a small distance from the origin, the interface collapses at the origin and the dynamic 
solution relaxes to Q2 in (37). In particular, |Q(r,t)p ^ |/i^ uniformly on the cylinder for large 
times. Figures [T^ and [14| show the time snapshots of the spatial distribution of eigenvalues, and the 
time evolution of the eigenvalues at r = 0, showing the convergence to a uniform uniaxial solution 
throughout the cross-section. Figure [TS] shows snapshots of the relaxation of one component of Q 
to the corresponding Q2-component. 

We also study how the initial non-planarity (as measured by e) affects the characteristic relaxation 
time to Q2. We observe that the modulus, |Q|^ (0,t) jumps abruptly from zero to at some 
critical time. Let be the first time for which 


and we associate with the loss of interface structure. Figure [TT] plots as a function of — logiQ e 
for various uq and vq^ and we find that oc — log^^gc. This can give quantitative estimates for 
the real-time persistence of isotropic-nematic interfaces and their experimental relevance for model 
problems with non-planar initial conditions. 


5 Biaxial boundary conditions on a disc 

The Dirichlet conditions in Sections and are uniaxial minima of the bulk potential, fs and are 
referred to as minimal boundary conditions. In this section, we employ a biaxial planar boundary 
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Figure 15: Qii(r,t) — ^(4x^/(l + — 1/3) on the cross section of the cylinder at 2; = 0 for 

initial condition ( [^ with uq = 0.6, vq = 0.4 at t = 0, t = 0.001 and t = 0.6. 



Figure 16: Interface position of dynamic solution (58) with initial condition (57), with uq = 0.6, 
Vq = 0.4 for L = 0.05 (grey) and L = 0.01 (black) compared to motion by mean curvature (blue). 
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Figure 17: Time (t*) at which the interface is lost for various e and uq and vq] uq = 0.1 (purple), 
Uq = 0.4 (blue), Uq — 0.5 (green), uq = 0.6 (orange) and uq = 0.9 (red), 'L’o = 1 “ '^o initial 
condition ( [^ and L = 0.05. 
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condition that is not a minimum of fs at the nematic-isotropic transition temperature, given by 


Q> = (n, 


ni — m 0 mj 


(63) 


This boundary condition is maximally biaxial with a zero eigenvalue i.e. trQ^ = 0. We refer to 
as a non-minimal Dirichlet condition. From [29], we expect that LdG energy minimizers, subject 
to a boundary condition of this form, converge in an appropriate sense, to a limiting harmonic map 
that is a minimum of the bulk potential almost everywhere and develop a boundary layer near 
r = 1 to match Q^, in the vanishing elastic constant limit. 

We study two dimensional (2D) and 3D dynamic solutions separately. A 2D solution is a symmetric 
and traceless 2x2 matrix m and in such cases, we study maps from the disc to a 2D target space, 
with just two degrees of freedom. A 3D solution is a symmetric, traceless 3x3 matrix and in such 
cases, we study maps from a 2D domain into a five-dimensional target space. 

We start this section with a discussion of the 2D case; 2D Q-matrices have trQ^ = 0 and the 
corresponding evolution law simplifies to 


Qt = LAQ-^Q-C|QpQ. 

The simplest 2D dynamic solution, consistent with ( |6^ , is 

Q(r, t) — s(r, t) (ni (g) ni — m (g) m) 


(64) 


(65) 


with + y^. It is simple to check that the gradient flow model (64) admits a solution of the 

form (65) if the function s{r,t) is a solution of 


ist = 

with fixed boundary conditions 


+ 


Sr 45 

s / B"^ 

r j 

L \27C 


B 


+ 2Cs^ 


( 66 ) 


s{0,t) = 0 s{l,t) = — 


(67) 


for all t > 0. The evolution equation ( |6^ is simply the gradient fiow model for the functional 

/ 


i[s] := r 

Jo 


2s^ 

dr J r 


r 

+ - 


L \27C 


5^ + 


( 68 ) 


and given a smooth solution, s{r^t) of (66)-(67) with suitable initial conditions, ([6^ is the unique 
2D solution. Further, the 2D potential has an isolated minimum at s = 0 (see (|68[)) and hence we 
expect that any dynamic solution of (64) has an outward-propagating interface that separates an 
almost isotropic core around r = 0 from the Dirichlet boundary condition at r = 1. The interface 
equilibrates near r = 1, followed by a sharp boundary layer to match the fixed boundary condition. 
Next, we present some heuristics for 3D dynamic solutions with a planar initial condition of the 
form 

Q(r, 0) = s (r, 0) (ni (g) ni — m (g) m) , (69) 

where 5(r, 0) has an interface structure i.e. 


5(r, 0) 



r <ro 
ro < r < 1 


(70) 
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Figure 18: Radial profile of Qii{r^t) — — 1/3) for 0 

ro = 0.5), at t = 0, t = 0.001 and t = 0.25. 


0 for initial condition (69) (with 


for some ro G (0,1). Based on the analysis in Section Em all dynamic solutions remain planar 
with Qis = Q 23 = 0. Hence, we expect Q(r, t) (for long times) to have an isotropic core at r = 0 and 
Q(r,t) converges to Qi (defined in ([37|)) everywhere away from r = 0, with a boundary layer near 
r = 1 to match the Dirichlet condition. However, we speculate that there is a second scenario for 
non-minimal boundary conditions as in ( [6^ which is not observed for minimal boundary conditions 
as in Section]^ If 1 — ro is sufficiently small, 3D solutions may exhibit an outward growing isotropic 
core since the isotropic phase is also a minimizer of the bulk potential and this scenario may be 
energetically favourable. All dynamic solutions with non-minimal boundary conditions develop a 
boundary layer near r = 1, which has an energetic cost. However, solutions with minimal boundary 
conditions (as in Sections and do not have boundary layers near r = 1 and in such cases, it is 
energetically preferable to either have a localized core of reduced order near r = 0 (as for planar 
3D solutions) or to have uniform order throughout the disc (as for non-planar 3D solutions). 

Thus, we expect sharp contrast in the behaviour of 2D and 3D solutions. The 2D dynamic solutions 
have little nematic order being largely isotropic or close to isotropic, except near r — 1 and 3D 
dynamic solutions are largely uniaxial with perfect nematic ordering (for at least a range of values of 
ro), except near r = 0 and r — 1. It is unknown whether 2D dynamic solutions could be physically 
relevant for severely confined cylinders and if so, what happens at the critical transition between 
2D and 3D dynamic solutions. In the next section, we present numerical results to validate our 
heuristics above. 


5.1 Numerical simulations with biaixial boundary conditions 


We numerically solve ( |27{|31| ) on a disc with the Dirichlet condition (63) on r = 1. 
parameter values are as in Section |4.2[ 


All other 


(69) with s(r) = ^ 


1 + tanh 


with L = 0.01. The planar initial condition is as in 
In the case of ro = 0.5, the solution quickly becomes 


almost uniaxial in the interior by developing an inwards-propagating well-defined isotropic-nematic 
interface. This is illustrated by a plot of eigenvalue evolution in Figure The solution converges 
to Qi(r,t) (see Equation([37|) away from r = 0 for long times, modulo the core of reduced order 
near r = 0 and a thin boundary layer near r = 1, as displayed in Figure [T8| This is as expected 
from the numerical results presented in Section 

Next, we consider ro = 0.92 and observe a different behaviour; the interface evolves so there is 
a thin boundary layer near r = 1 with a large almost isotropic core in the interior. This is best 
illustrated with radial profiles of |Qp as seen in Figure 20, 
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Figure 19: Eigenvalues of Q(r,t) for initial condition (69) (with tq 
t = 0.25. 


0.5), at t = 0, t = 0.001 and 
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Figure 20: Radial profile of |Q(r, t) f for initial condition (69) with ro = 0.92, for t = 0 and t — 0.25. 


We compare the 3D solutions above with 2D solutions for the same system ( [6^ . We work with 
planar initial conditions ( [6^ , with 5(r) as before. The initial interface grows rapidly to yield an 
almost entirely isotropic interior with a thin boundary layer near r — 1. This is illustrated by the 
eigenvalue evolution in Figure thus corroborating our heuristics and analytical reasoning in the 
previous section. 

6 Conclusions 

We focus on the gradient flow model for the LdG energy on prototype geometries, such as a droplet 
and a disc, with Dirichlet boundary conditions and various initial conditions at the nematic-isotropic 
transition temperature. In Section we consider the model problem of a 3D droplet of radius i?, 
with radial boundary conditions. In the case of uniaxial radially symmetric initial conditions with 





Figure 21: Eigenvalues of Q(r,t) in the 2D model, for initial condition (69), at t = 0, t = 2 x 10 ^ 
and t — 2 X 10“^. 
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an isotropic-nematic front structure, we adapt Ginzburg-Landau methods from main] to prove 
that the isotropic-nematic interface propagates according to mean curvature in the ^ 0 

limit. However, the qualitative dynamics seem universal for a large class of radially symmetric and 
non-symmetric uniaxial and biaxial initial conditions and the long-time dynamics are determined 
by the classical RH solution, which has been numerically demonstrated to be a global LdG energy 
minimizer in this regime. 

In Sections and we focus on dynamic solutions on a disc. Our results are largely numerical 
and complemented by heuristics and analytical reasoning. We demonstrate how a choice of planar 
or non planar initial condition can influence the long-time dynamic behaviour. Planar initial 
conditions generate planar dynamic solutions with an isotropic core at the centre for all times 
whereas non-planar solutions follow the planar dynamics for a length of time, before relaxing into 
an uniaxial state of perfect order for long times. In Section we look at non-minimal boundary 
conditions. Non-minimal boundary conditions allow for dynamic scenarios outside the scope of 
minimal boundary conditions and since minimal boundary conditions are an idealization, non- 
minimal Dirichlet conditions can be physically relevant too. 

The long-time dynamics can be understood in terms of local and global minimizers, or in some 
cases critical points, of the LdG energy. In cases where the LdG critical points exhibit an isotropic- 
nematic interface, this interface may be localized with little effect on global properties. Our nu¬ 
merical results show that a large class of physically relevant LG model problems can exhibit a 
well-defined isotropic-nematic interface for a length of time (see Fignres and [T^ and these results 
give insight into how boundary and initial conditions can be used to yield either largely disor¬ 
dered or ordered nematic profiles. A natural next step is to rigorously analyze front formation and 
propagation with generic non-minimal boundary conditions and with more general LdG energy 
functionals, including those with a sixth order bulk potential that allow for biaxial minima. We 
will report on these developments in the future. 
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